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A theory of self-assembly of surfactant molecules into micelles and bilayers is critically examined and 
extended to include vesicles and microemulsions. The notion of hydrophilio-lipophilic balance is 
quantified. The theory gives a unified account of type, size and shape of the aggregates which form under 
various conditions. Observed trends due to change in salt concentration, temperature and oil type, and 
due to the addition of cosurfactants, are correlated and emerge from a simple global framework. 

1. INTRODUCTION 

Until a few years ago the possibility that all observations on association colloid 
systems could ultimately be handled by a single theoretical framework seemed 
remote. It became less so following attempts 1,2 to extend the ideas of Tanford 3 and 
others 4 on dilute micellar aggregates to larger surfactant associations like cylindrical 
micelles, vesicles and bilayers. The main point of departure lay in quantifying the 
part played by molecular geometry (packing) in determining allowed structures. It 
was an old idea 5 which had been allowed to lie fallow. And it worked. Theory does 
appear to be on the right track. While there are gaps, parts of the jigsaw puzzle have 
been filled in more or less satisfactorily for dilute surfactant solutions. Certainly 
many of the physical properties of micelles and vesicles like size, shape, c.m.c. and 
polydispersity appear to be accessible without a detailed knowledge of the complex 
intermolecular forces involved. 

Our purpose here is twofold: (1) To attempt to define better and explore some of 
the basic assumptions which underlie ideas presently extant. (2) To see how these 
ideas might be extended to include multicomponent systems (microemulsions). 

From a pragmatic point of view, one main aim of studies in the subject must surely 
be: to elucidate the phase diagrams of water-surfactant (and cosurfactant)-hydro- 
carbon mixtures; in particular to identify which structures form, when and why; and 
as a corollary: how to maximise solubilisation of oil in water, or water in oil, with a 
minimum surfactant (cosurfactant) concentration. 

This aim is ambitious, and the problem of such complexity that, ttfparaphrase and 
borrow a remark made by Stillinger: 6 "it is essential to maintain a respectable 
balance between the sterile intricacy of formal theory and the seductive simplicity of 
poetic 'explanation"'. Before beginning our study it may be useful to expand this 
dictum. In attempting to make a theory there are two extreme approaches. A 
fundamental treatment using statistical mechanics which takes into account complex 
surfactant molecule interactions in water is possible in principle. However, even the 
hydrophobic interaction between two small molecules in water is still a matter of 
dispute/ Further, the simplest prototype for aggregation, the probem of nucleation 
(and consequent phase transition) in a van der Waals gas, is an open subject. 
Moreover, the high road via statistical mechanics is necessarily so complicated that 
physical insight tends to be wholly obscured, 
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At the other end of the spectrum one can avoid detailed models as far as possible 
and search for a unified thermodynamic picture like that of Tanford, 3 and of 
Israelachvili et al l for micelle and bilayer formation. But thermodynamics is 
tautological and can go so far. At a certain point some details of molecular 
interaction must be invoked. Our guiding principle in attempting to steer a middle 
course between these two extremes will be that these details must be minimal. 
Otherwise, with too many (unknown) parameters, theory tends to become aji 
exercise in curve fitting, a numerical game which loses predictive capacity and 
credibility. 

An immediate consequence is that language must be used with some care. The 
problem is here doubly compounded and confounded by the fact that words like 
amphiphile, hydrophobic, hydrophilic, lipophilic, aiggregate, micelle, are either 
. anthropomorphic in origin (disguising their complexity) or intuitive, ill-defined, and 
are so familiar that we tend not to question their meaning. This question of 
definitions will plague us throughout and is unavoidable. A result is that our essay 
has elements of schizophrenia, with necessary appeal to formal statistical mechanics 
(relegated to an appendix) interspersed with (hopefully) an occasional insight 
gleaned from intuition and simple models. 

2. THERMODYNAMICS OF DILUTE SOLUTIONS OF 

SURFACTANTS 

Dilute micellar systems have been the subject of an extensive literature which we 
do not wish to recapitulate in detail. Given certain "reasonable" assumptions they 
are broadly speaking understood. However, conflicting opinion on questions like the 
stability of vesicles and the delicacy of experiments on microemulsions, to name two 
areas of many, will demand re-examination of foundations as we proceed. 

For the moment it suffices to assert that a dilute solution of surfactant molecules 
can be considered to consist of solvent plus monomers, dimers, trimers, . . . and larger 
allowed aggregates (micelles, vesicles, liposomes, . . .). The concentration is assumed 
to be so low that aggregates can be considered to be non-interacting. The dis- 
tribution of aggregates is then determined by the law of mass action, eqn (1), 

>*+|^ (1) 

where the chemical potential of an aggregate of size N has been written as 
NfiM+kT In (Xn/N) and Xs is the concentration (mole fraction, volume fraction, 
, . . .) of surfactant molecules in the JV-aggregate. The theory also allows for the 
possible formation of infinite aggregates, i.e. separate phases, a problem which we 
defer. The glib assertion eqn (1) represents a beginning to a chemist and an essential 
stumbling block to a physicist who can go no further without questioning origins. If 
aggregates of a given N were distinct, identical well-defined chemical species there 
would be no problem, apart from the vexed question of concentration units. They 
are not: even within a given TV-aggregate, if such can be defined, there exist an 
infinite diversity of conceivable shapes or configurations which the association of 
surfactant molecules could take up. Implicit in eqn (1) is the understanding that for 
any N a shape of minimum energy exists and is overwhelmingly more probable than 
its fellows. These problems are addressed in the Appendix. If the arguments there 
advanced are accepted we return to eqn (i) which can be rewritten in the form 

&=X?e^ (2) 
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All aggregates can occur at any concentration, albeit with infinitesimal probability 
even below the critical micelle concentration (cm.c.). Above the c.ra.c, defined by 
X\ = l£>i Xm Xx increases slowly with concentration. It is shown in ref . (1) that if 
ixifis sharply distributed about some N then the distribution of aggregates peaks at a 
value of N just less than the TV with minimum ft* and is also sharply distributed. 
Otherwise pronounced polydisperity may occur {e.g. for long, cylindrical micelles). 
Thus reduced to bare bones, it can be seen that the use of the word "theory*' is 
dubious. We have simply characterized the observation of micelles, and have shown 
via the Appendix that the law of mass action is an appropriate vehicle for this 
characterisation. 

The entropic term (kT/N) In (X N /N) in eqn (1) has considerable nuisance value. 
For that reason it will usually be helpful in dealing with finite aggregates to adopt a 
simpler basis for our subsequent development This is the so-called pseudophase 
approximation wherein one drops the entropic term. In the pseudophase approxi- 
mation no micelles occur below the cm.c. (This is now the value X\ot monomer 
concentration for which v>N = li>\ + kT\n Xx) Abbve this cm.c. all additional 
surfactant molecules form micelles or whatever aggregate has the minimum fi° N . No 
other aggregates form until activity coefficients, Le. interactions between aggregates, 
become significant. Keeping in mind possible complications due to phase transitions 
and interactions the strategy is then to compare the chemical potentials of different 
aggregates to see which has the minimum ^. 

.. . THE FORM OF fls 

• In essence the simplest version of current theory 1 assumes a form 

^N = MN+MN^-M^+ t< packingterm ,, (3) 

= "bulk term" + "surface term" + "curvature term" + "packing". 
Its decomposition is as follows (c/. fig. 1). 




Fig. 1 .—-Schematic representation of model spherical micelle of radius R> 

Bulk Term : This is a constant term, the same for all aggregates, which measures 
the hydrophobic free energy of removing hydrocarbort tails from water into an 
assumed oil-like phase made up of all tails which form the micelle interior. The 
interior is assumed to be fluid in estimating this free energy transfer. 

Surface Term : This includes a termt ya to allow for the fact that hydrophobic tails 
still have some contact with water, where a is the area per surfactant molecule, and y 
is an interfacial tension. 



f It is a trivial decoration to exclude that fraction of the "oiT-water interface taken up by head groups. 
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Opposing this attractive energy is a term to account for repulsive head-group 
interactions. These interactions can be due to steric repulsion, hydration, electro- 
static and other forces. They are not yet quantified. If electrostatic in origin the 
simplest phenomenological description would give a term oc 1/a. The total surface 
contribution would then be written 



which takes its minimum value 2yao at an optimal area ao per head group. This form 
for the repulsive contributions to fi% cannot be taken literally, and its precise form 
is immaterial. Any mathematical description of the surface contributions which 
recognises that several competing forces will lead to an optimal area a 0 per head 
group will lead to the same conclusions. 

Curvature Term : If as an aid to analysing data on ionic micellar solutions we were 
to persist with an electrostatic model for head-group repulsions, 1 then for a curved, 
rather than planar, surface the term a in yal/a would be replaced by a(l+D/R) 
where R is the mean radius of curvature and D the Debye length. An alternative 
way of visualising a mechanism for curvature effects is to imagine that the centre of 
action for head-group repulsion is displaced a distance. I) from the "oil'-water 
interface. Then a would be replaced by a (1 +D/R) 2 .. We choose the first form again 
only to have some plausible mathematical realisation of the undoubted existence of 
curvature energy. 

Packing Term : The assumption that the interior of an aggregate is fluid-like and is 
to a first approximation incompressible has an immediate consequence, provided we 
admit that aggregates can obtain no holes. [The occurrence of an interior, vacuum of 
water-filled region inside the (oil-like) interior of an aggregate would result in a large 
unfavourable increase in free energy, which possibility is excluded from 
consideration]^ This can be taken into account if we assume fih-^ hence a 
packing criterion is violated. For spherical and cylindrical micelles this criterion is 
R < l c where / c is a critical tail-length which in ref . (1) was taken to be ca. 80-90% of 
the fully extended chain length for bilayers. The packing criterion is clearly an 
extreme, albeit useful, oversimplification. We shall see later that it can, and indeed 
must be relaxed. The melding of the two notions, of a fluid-like interior for the 
micelle and of "packing", is at first sight contradictory. However, the two notions 
can be shown to be compatible in a first-order theory. 8 

In addition to the contributions above to p* N there will be others with increasing 
concentration due to interactions. 

3. GEOMETRICAL CONSIDERATIONS AND ZEROTH-ORDER 

THEORY 

Possible candidates for aggregates can now be examined. For surfactant-water 
systems these are spherical micelles, non-spherical micelles (globular and cylindri- 
cal), vesicles, liposomes, bilayers, and for oil-surfactant-water systems spherical 
drops, normal or inverted (water dispersed or' oil dispersed, respectively). We 
assume that the surfactant molecule can be divided unambiguously into a head group 
and tail and that the partial volume v of the tail is prescribed so that for an aggregate 
of number N the tail region volume is Nv. 

The location of the interface which defines the radius R is a difficult question with 
different conventions adopted by different authors. 1,9 However if the theory is not to 

t Note added In proof 9 . Experimental observations which have sometimes been Interpreted as implying 
water penetration in the hydrophobic core of micelles are completely explicable in terms of the sampling 
of the surface groups in the chains (D, W, R, Gruen, to be published). 
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depend on specific details, the convention must also n t be critical. Take the 
interface to be the surface which bounds the volume occupied by hydrocarbon tails. 
Then for spherical micelles (fig. 1) we have iirR 3 = Nv and a head-group area defined 
by 4irR 2 =Na. The packing criterion imposes the condition that the average chain 
length is R ^ / c . 

By zeroth-order theory we mean one in which a - a 0y the optimum area for a 
bilayer, independent of the size and shape of the micelle. (Head-group curvature is 
assumed to be a small perturbation.) Because the zeroth-order theory prescribes 
a = a 0 , the dimensions of the spherical micelle are completely determined (R~ 
3v/a 0 ) r i.e. R>l c . Then if v/a 0 l c >h "the optimum-energy spherical micelle" is 
disallowed by packing because then R > / c . On the other hand, for R ^ spherical 
micelles allowed by packing are energetically disfavoured (because a > ao). In 
higher-order theories a will be determined by optimisation of fi# and the optimum 
spherical micelle (R ^ / c ) will be compared with other shaped aggregates. 

Problems of packing are extremely complex for shapes of low symmetry. For 
shapes of high symmetry the packing criterion is easy to visualise. In other cases the 
simplest condition is that no interior point can be less than a distance / c from the 
surface. This criterion leads to the conclusion that only oblate shapes should 
occur. 1,10 In practice this is not so, and large prolate or rod-shaped micelles are 
observed. 1,28 

To overcome this difficulty in ref. (1) the authors imposed a plausible restricted 
(local) packing criterion. As there shown, the pseudophase approximation breaks 
down, and large polydisperse cylindrical micelles form. Their size distribution 
depends on concentration. So as not to distract from our main theme it will be 
sufficient to treat non-spherical micelles as if they were infinite cylinders. For 
cylinders R = 2v/a 0 so that packing requires v/a 0 l c <2- 

So far for v/a 0 I c < h the zeroth-order theory admits either a sphere or a cylinder. 
To determine which is the favoured aggregate we can have recourse to two 
arguments. In ref. (1) appeal was made to translation entropy which favours the 
smallest aggregate, Le. a sphere. However, we can invoke the head-group curvature 
energy which also prefers the sphere because of its greater curvature. 

Note that micelles are forbidden for v/a Q l c >2. In this region we have to search 
for different allowed structures. These could be, e.g., vesicles, liposomes or multi- 
lamellar aggregates. For the moment consider only vesicles and single bilayers* In 
discriminating between these, the head-group curvature again favours the smallest 
aggregates as does entropy, le. the smallest vesicle allowed by packing is preferred. 

The zeroth-order theory also allows vesicles for v/aol c <\ However, it will be 
shown in section 4 that the preferred vesicles have physically absurd internal 
parameters which necessitate the imposition of a further packing constraint. Once 
this is taken into account cylinders become the energetically favoured structures for 
l<v/a 0 lc<i 

With this proviso, the overall prediction is that we obtain 
spherical micelles v/a 0 l c <3 

cylindrical micelles k<v/a Q l c <2 (5) 
bilayers (or vesicles) 2<tvW c <l 
inverted structures v/aol e >l. 

Apart from the question of internal packing alluded to above, there are other 
problems with this theory which require detailed examination: We have no reason to 
assume that curvature effects due to head-group interactions are necessarily small, in 
the sense of a perturbation. If they become large, the zeroth-order theory breaks 
down in any event because the head^group area may no longer be fixed to first 
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approximation as equal to that for a (planar) bilayer. We have no reason to assume 
that curvature is predominantly due to head-group interactions. There may well be 
curvature effects due to hydrocarbon tails of comparable magnitude which could 
even take a different sign (vesicle-forming surf actants generally have relatively large 
hydrocarbon volumes). " 

In addition the experimental situation is still unresolved. Some authors 
believe that vesicles are thermodynamically unstable (the stable state being the 
lamellar phase). Monodisperse vesicles of size given by the zeroth-order theory 2 can 
be produced through prolonged sonication. This is generally taken to indicate that 
energy is required for their production. When monodisperse vesicles are so 
produced and cooled below the gel temperature, they become unstable and grow in 
size, eventually forming lamellar particles. 13 On reheating, whatever size dis- 
tribution existed at the lower temperature is apparently "frozen" and seems to be 
kinetically stable. 13,14 Again it is claimed that with other preparation techniques any 
size vesicle can be made to order. 15 Not all these different distributions of vesicles 
can be the thermodynamically stable state. All that can be inferred is that vesicle 
dispersions are very slow to equilibrate, so that it is difficult to ascertain the 
equilibrium situation. 

On the other hand, on addition of excess water, the laitiellar phase disperses 
spontaneously to form multilamellar particles or liposomes 16 ' 17 (e.g. lecithin) or 
vesicles 1748 {e.g. phosphatidyl serine). This would seem to suggest that the lamellar 
phase may not be stable in a dilute solution. 

Another problem which we have not yet addressed is that of interactions between 
bilayers. The question of stability of vesicles cannot be resolved theoretically 
without proper account of interactions. For all these reasons the problem of vesicles 
needs to be reconsidered with care. ' 

4. VESICLES 

4.1 PACKING IN VESICLES 

We consider first consequences of packing conditions which apply under the 
assumptions of (1) constant volume v per surfactant molecule and (2) constant head- 
group area ao for a surfactant molecule. In general a vesicle {cf. fig. 2) of aggregation 
number N is described by four parameters of the set v, a e , a i9 U> U, R*> Ri> W N ei N lf 
where the subscripts e, i refer to external and internal parameters. Of these we take 
the first four v, a c -a t - a 0> U to be given. We have 

y(*e-*^NeV 4irRl=Nid 0 , R t -R m =le (6) 




Fio, 2.— Geometry of a model vesicle, 
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whence 

U 3 l3v 3 1 v 4 

jR e 2 v a 0 l t 4 3 aol c 

As a/a 0 fe increases from 5 to 1, RJU increases from 1 to 00, and R c increases from 
3i?/a 0 tooo. Hence for fixed v,a 0 the smallest vesicle is given by the largest f c . Thus, 
since l c ** / c , the smallest vesicle corresponds to l c - l c (provided v/a 0 lc>3)- 

III the following we determine the dimensions of the vesicles as a function of the 
simple parameter v/a 0 l e . The numbers listed in table 1 follow from the relations eqn 
(6) and (7) and the further conditions 

y(R* T *?) = N* 4ir{Rl+Rf)=Na 0 (8) 

whence if x = RJRe, y = 3v/(a 0 R e ) 

Jjp=y, Q<x,y<l. (9) 

Formally the model allows vesicles (with / e = / c ) for v/aok > 3 . However, as can 
be seen from table 1 , in practice for 5 < v/a 0 l c < 2 these vesicles are excluded. In this 
region the number of surfactant molecules "packed" on the inside of an optimum 
vesicle is less than a whole number, which is inanfe. Clearly some additional packing 
criterion must be invoked since the values of U for the mean length of internal chains 
are ridiculously small for v/aoU<h Within the framework of our zeroth-order 
theory (fixed head-group area) this condition would have the form h**al Ci 0<a < 1 
where a is a parameter ultimately related to details of chain-chain interactions. At 
the very worst we must have /j>2-3 A say, roughly the length of a CH 2 bond. 
Assume that such a condition is invoked. Then the smallest vesicle allowed by 
packing in the range \<v) Wc<? will correspond to some / e , 5<v/a 0 U<l. 
Depending on internal packing constraints, ie. the value of a, such vesicles may have 
small aggregation numbers. Thus if entropy was the sole deciding factor, vesicles 
would always be favoured over the large cylindrical micelles observed in practice. 
[As shown in ref. (1) cylindrical micelles once formed tend to grow without limit.] 
Therefore another factor must decide the issue. We now show that once curvature 
energy, however small, is taken into account cylinders are almost certainly always 
energetically favoured oyer vesicles in the range |< v/a 0 l c ^h 

4.2 CYLINDERS VS. VESICLES 

Plausible simple forms which take account of curvature energies due to head- 
group repulsion have been alluded to in section 2. These will be discussed further 
below (section 4.3) and are clearly model-dependent. Thus an electrostatic model 
for the surface contribution to fi %r [eqn (3)] would replace the term yal/a for a planar 
aggregate by yal/a(l+Di/R) for a spherical aggregate. A different model which 
recognised head-group size would replace this term by yao/a(l+D2/R)*. 
However, provided curvature represents a weak perturbation, linearisation is 
allowed. All models then lead to the same approximate form y(al/a)[l-D/R] 
where 1/R = £[(l/#i) + (l/l?2)], #1 and J? 2 being the principal radii of curvature of 
the aggregate, and D is some unknown constant characteristic of a given surfactant. 
If curvature due to head-group interactions dominates over chains, which is implicit 
in our zeroth-order theory, D must be positive for a convex surface and negative for a 
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concave surface. We now compare the energies of cylinders and vesicles (with the 
same v and a 0 ). For cylinders we have 

where R c is radius of the cylinder. By contrast, for vesicles we have 

/tt N (vesicle) = (11) 

N 

where 

Hence /* s (vesicle) can be written 

S/ w _ Da 0 /N c NA _ 3Pt> 

where the last form follows from elementary geometry. Comparison of eqn (10) and 
(13) shows that cylinders will be energetically favoured over vesicles if and only if 

K* + W*f»(?H) 2 ; or §f±&£ (14) 

We now establish conditions under which the inequality holds. To do this we 
write r c « RJ U> r% » Ui// C> x = RJR { so that the inequality (14) becomes 

r«+rf 4t> 1 + x 2 4i> 

- — L >— r or - 5 >— r- ( 15 > 

Now r c , ri and x are all functions of v/a 0 l c - In fact from eqn (7) 

r e \2 » a 0 / e 4/ 

and x is determined by eqn (9). Hence the condition eqn (15) can also be written as 

-^^>i (17) 



(t-x)a-x 3 ) 3 

The solution of 

(1-*)(1-* 3 ) 3 



2\2 



a+x<) 

then determines vesicle parameters for which the inequality is satisfied. The 
required root is x = 0.1896, corresponding to which value we have from eqn (9) and 
(16) t>/flo/e = 0.4883, t>/a 0 /i = 2.104, J* e //e = 1.528, Ri/U = 0.289, /,// o =0.232. 
Therefore, vesicles which satisfy the conditions 

-^r> 0.4883 or equivalehtly -*V< 2.104 or 0,2321 (18) 

Clo'o a o'l 'e 

will be energetically less favoured than the cylinder. 
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We now. show that all vesicles satisfy this condition given the minimal packing 
condition of section 4. 1 . If v/a 0 l c > 0.4883 the condition is satisfied since / c * 1* This 
is not a particularly interesting result since we know that cylinders are forbidden by 
packing in any case for v/a o l c >0.5, On the other hand if v/a 0 l c < 0.4883, the 
external packing condition U^h is not sufficient to discriminate, and a packing 
condition on internal chains [cf. eqn (18)] is required. Let us impose the requirement 
jfi// tt >0.2321. Then certainly ljl c > 0.2321. Thus for v/a 0 l c < 0.4883, the condition 
h/U> 0.2321 is sufficient to ensure that cylinders are again the favoured structure. 
Typically for a lipid of l c « 20 A, the mean internal chain length of a vesicle which has 
an energy equal to the allowed cylindrical structure would be as low as 4,6 A. This is. 
improbable. The corresponding vesicle would be an absurd structure, with aggre- 
gation number 26, N t « 2. 

To sum up: Within a framework which admits the zeroth-order theory as a first 
approximation, i.e. provided curvature is a perturbation, cylinders are always the 
favoured structure for 3 < v/ aol c < \. The imposition of a weak packing criterion on 
inside chains is necessary. This criterion cannot be spelt out precisely without a 
theory of chain interactions presently not available. However, the very strong 
requirement l x > 0.232/ c appears to be eminently reasonable. Entropy which favours 
the smaller vesicles over cylinders has been ruled out as a major factor and curvature 
always dominates. {To quantify this conclusion roughly compare the entropic 
contribution to fix of eqn (1) with curvature energies. Take X N = 1 as an extreme 
upper bound. Then even for a small vesicle 500, R e **.50 A, y 50, 

if 

For D as small as 1 A, the smallest possible characteristic length which has much 
meaning, the curvature contribution to would be [cf. eqn (13)] >4Dv/R e =* 
20xl0" 2 ikr.} 

5. FURTHER CONSIDERATIONS ON VESICLES— 
NON-LINEAR EXTENSIONS 

5.1 VESICLES VS. BILAYERS 

The zeroth-order theory of section 4.1 suggests that vesicles would be the 
preferred structure in the range §<t>/ap/ c <l. However, we have shown that the 
weak packing criterion for hydrocarbon chains / e < / c is not sufficient. Recognition of 
the role of an internal packing constraint can dramatically alter the picture. We first 
investigate vesicle parameters in the zeroth-order approximation with the constraint 
h>al c . The strong condition / s 2s 0.25 / c established in section 4.2 provides one 
bound on a. Evidence from multilayer studies suggests that the upper bound on a, 
a *z 0.5, is a reasonable guess. 19 Vesicle parameters for a = 0.5 are listed in table 2. 
These results were obtained numerically by solution of eqn (6)-(9) with the addi- 
tional constraint U >al c > 

Note that for a = 0.5, as v/a Q l c decreases from unity (bilayer formation favoured) 
vesicle size diminishes to around v/ ado- 0.7. Thereafter the optimal vesicle size 
increases until at v/a 0 l c = ot = 0.5 bilayers are again the favoured structure, The 
turning point can be read off from table 1. Thus if 

7i R m -R\ - _ 
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Table 2.— Effect of internal packing constraints on vesicles when 

a = a 0 ,a=0.5, u = 1063A 3 , / C =20A 
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. 10 


560 


130 


(0.65) 


82 


61 


33 


18 


10 


580 


170 


(0.59) 


90 


79 


54 


14.3 


10 


880 


420 


(0.56) 


94.5 


98 


76 


13 


10 


1300 


770 


(0.54) 


98 


137 


115 


11.9 


10 


2430 


1720 


(0.50) 


106 


00 


CO 


10 


10 


00 


oo 



Vesicle parameters with v = 1063 A 3 (corresponding to egg lecithin 1 ) and U taken to be 
20 A, with varying head-group area a 0 . ot is taken to be 0.5. 

corresponding to rather "stiff" chains the turning point is v/a Q l c =0&. Bilayers 
would be formed for 0.5 <t;/a 0 / c < 0.65 and revert to vesicles again for v/aoh> 
6.65. This remarkable phenomenon is purely a consequence of geometry and is not 
universal. Thus for very fluid chains, say a =0.25, vesicle size decreases mono- 
tonically until v/a 0 l c = 0.5, below which value cylinders are the favoured structure. 

5.2 NON-LINEAR EXTENSIONS 

Subject to those complications just discussed which arise due to internal packing 
constraints, both entropy and head-group curvature tend to favour vesicles as 
compared with bilayers. We wish to explore how that conclusion holds within a 
wider framework. That is, we now admit the possibility that curvature effects are 
large, so that the assumption a = a 0 of the zeroth-order theory would no longer be 
tenable. A theory which covers this situation is necessarily model dependent and 
non-linear, depending on unknown details of molecular interaction. 

Perhaps the simplest mathematical description of head-group interactions for 
non-ionic surfactants is via the capacitor model: Imagine that D is the average 
distance between positive and negative charges which form a dipolar head group {i.e. 
qD is the dipole moment where q represents unit charge). Then for a bilayer we 
should have 



s , WP l27rg 2 D /im 



where e is the "capacitor" dielectric constant. (A similar model might be used to 
describe ionic surfactants with D the Debye length.) Corresponding to eqn (11), we 
now have for vesicles 

. (20a) 
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where 

N t ^(Rl -*»); W-^(Ri-Hf); N=N e +N-,. (21) 
3d 3t> 

Recall that J? c , -R m and i?i are external, middle and internal radii, respectively, with 
AT e and iVi the number of external and internal surfactant molecules- Obviously 

Rc>Rm>Ri>D. 

Optimisation of fi y C$ is a tedious but straightforward numerical problem. Typical 
results are given in table 3, which should be compared with table 2. (For definiteness 
we take 4 = 4.8xl0~ lo e.s.u., *=78, y^SOdyncnT 1 so that ab = v(372U)A , 
where D is measured in A.) The particular algorithm used was that developed by 
Brent. 20 



Table 3.— Typical vesicle parameters in capacitance model for head- 
group interaction with a = 0.5, v = 1063 A 3 , / c - 20 A 



D/A 3 v/aok 


flo/A 2 


ajk 1 


a./A 4 Rjk % 


K./A 


h/k 


/,/A 




Ni/A 


7.55 1 


53 • 


S3 


53 


CO 


oo 


20 


20 


00 


00 


8 (0.97) 


54.5 


54.8 


55.3 


646 


607 


20 


18.6 


95 000 


83 000 


10 (0.87) 


61 


62 


68 


130 : 


96 


20 


13.7 


3430 


1730 


12 (0.79) 


67 


68 


85 


84 


53 


20 


10.4 


1300 . 


430 


14 (0.74) 


72 


68.4 


88.4 


83 


53 


20 


10 


1250 


390 


16 (0.69) 


77 


67.6 


89 


86 


55 


20 


10 


1365 


440 


18 (0.65) 


82 


68.4 


91 


92.5 


63 


19.3 


10 


1570 


550 


20 (0.62) 


86 


72 


94 


112 


84.6 


17 


10 


3150 


950 


.22 (0.59) 


90 


78 


98 


143 


118 


15 


10 


3 270 


1780 


24 (0.56) 


94 


84 


100 


200 


173 


13.6 


10 


5 800 


3 750 


26 (0.54) . 


98 


90 


102 


303 


281 


12.2 


10 


12 800 


9 680 


30 (0.50) 


106 


106 


106 


00 


00 


10 


10 


oo" 


OO 



An internal packing condition / } ^ al c is again necessary to obtain sensible results. 
Variation of a*, a\ from ao is considerable at larger head-group areas. Similarly, 
predicted vesicle sizes can vary by tip to a factor of 2 from the zeroth-order theory in 
extreme cases* Similar results emerge for different values of of. Despite these 
quantitative differences it does appear that non-linearity does not affect the main 
conclusions of the zeroth-order theory. 

We have to emphasize that whatever the nature of the forces between head 
groups (and chains) which give rise to curvature, these may well inject non-linear 
effects and precise quantitative predictions are not feasible at the present time. 
However, the general pattern which emerges from the zeroth-order theory seems to 
be on fairly safe ground. Moreover while the theory strictly depends on four 
parameters (v, ao, /« ot) the substantial predictions rely only on the surfactant 
parameter v/aoh and are independent of the detailed form of the forces. The vexed 
and delicate question of the stability of vesicles must remain open in general, 
depending as it does on chain stiffness (subsumed in the parameter a) and the 
competition between curvature and interaction between bilayers. The answer will 
differ from surfactant to surfactant and depends on temperature. 
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6. MICROEMULSIONS 

6.1 THERMODYNAMICS OF MICROEMULSIONS 

We now broaden our enquiry to include oil-water-surfactant systems. In the 
absence of surfactant, oil and water will separate into two phases (unless the volume 
fraction of oil is extremely low or very close to unity). In this case the oil aggregate is 
arbitrarily large, its size being determined by the absolute amount of oil present 
When surfactant is added, some surfactant molecules will go to the oil-water 
interface and lower the interfacial tension. With the addition of still more surfactant, 
several possibilities arise. Aggregates may appear in either the water phase or the oil 
phase. These aggregates may be micelles containing a little oil, dispersed in water or 
inverted micelles containing a little water, dispersed in oil; or large drops of oil (or 
water) with surfactant at the interface, dispersed in water (or oil). 21 * 22 It is possible to 
solubilize all the oil (or water). Another possibility is the occurrence of a three phase 
system: water, oil and so-called surfactant phase which contains all three 
components. 21 Once aggregates form, additional surfactant will be consumed in the 
formation of aggregates and will have little effect on the tension of the oil-water 
interface. We here use the term microemulsion to mean a thermodynamically stable 
dispersion of large oil and surfactant or water and surf actant aggregates. We make 
no distinction between large swollen micelles and microemulsions. Microemulsions 
are often monodisperse 23 and the size of the drops is determined by the relative 
proportions of oil, water and surfactant, not by the absolute amount of oil, in the 
solution. These are by definition one phase systems. 

The thermodynamic approach developed in the previous sections can be exten- 
ded to deal with these systems. Microemulsion drops can be treated in a similar way 
to pure surfactant aggregates and the distribution determined by the law of mass 
action. As for surfactants where a simpler picture is provided by the less rigorous 
pseudo-phase approximation; a simpler picture of microemulsions is also obtained if 
we ignore translational entropy of the aggregates. 

Recall that for a surfactant solution the aggregates which form (in the pseudo- 
phase approximation) are those with min /jln and the concentration at which they 
form (the cm.c.) is given by 

tih = V<i +fcrin*i=At. 

This result was derived as an approximation to the law of mass action, but it can be 
seen more directly as follows. If some of the molecules were arranged in non-optimal 
aggregates they could lower their total free energy by rearranging into optimal 
aggregates. If /x < fi the total free energy would be lowered by dissolution of some 
or all ,of the aggregates. We note also for optimal aggregates (defined by 
dAWaN = 0)that 

so that when (i = min fi^ the optimal aggregates are in chemical equilibrium with the 
surf actant monomers. 

In the case of multicomponent systems a similar argument can be invoiced. The 
determination of which mixed aggregates form is a little more complicated. Here the 
composition of the aggregates will be dependent on the overall composition of the 
mixture. (This is abundantly clear if almost all the solute molecules are aggregated,) 
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K G MtN ^(M+N)/i 0 MtN denotes the free energyt of an aggregate containing M 
molecules of oili and N molecules of surfactant then clearly out of all possible 
aggregates with the given composition X = N/(M+N) those which form must have 
the min filf tN . The minimization is taken with respect to shape and size, etc., at fixed 
composition. 

The composition of aggregates will be uniform if g = min filter is a concave 
function of X\ i.e. d 2 g/dx 2 ^ 0. However, if g is not a concave function the solution 
may contain a variety of aggregates of differing composition! (c/. fig. 3). 




X A 1 
X-Nt[M+N) 




0 Xj x A x 2 1 



Ms 



Fig. 3. — Diagram illustrating determination of aggregate composition: {a) g is everywhere a concave 
function of X. If chemical potential of one component is given (e.g. j*o)» composition X\dl aggregates 
and chemical potential of other component is determined, (b) g is not everywhere concave. For a 
particular value of /*o an d Ms the tangent touches the curve at two points. In this case a mixture of 
aggregates of compositions X u X 2 will occur, A system comprising aggregates of intermediate composi- 
tion would have a higher, free energy gA than the mixture g B . 

If we wish to determine the composition of the aggregates given the overall 
composition of the solution it is necessary to determine the concentrations, or 
equivalently the chemical potentials, of the monomers. The general problem is 
difficult. The easiest way to proceed is via an indirect route. We assume a value for 
the chemical potential of one species. This allows us to determine the corresponding 
value of the chemical potential of the other species. The composition of the 
corresponding aggregate and overall composition of the corresponding solution 
follow. The procedure is iterated until the required overall composition is achieved. 
One special case, indeed one of much interest, can be handled with relative ease: For 
a microemulsion solution in the presence of excess oil the chemical potential of one 
species, oil, is in fact prescribed. 

Suppose then that we are given the chemical potential /Lt 0 of the oil (i.e. of the 
monomers; and also of bulk oil when excess oil is present). We wish to show that the 
aggregates (and there may be more than one type) which can occur for this fio are 
those for which the expression fis(M 9 N) = {G MtN -M/jlo)/N \& a minimum. The 



t It is argued in the Appendix that (M + N)ftfcw is the Helmholtz free energy + P Q V flgg where P 0 is the 
(external) pressure of the solution. 

% The argument holds for water instead of oil. 

|| The argument here parallels the argument that the Gibbs free energy of a mixture of two components 
divided by the total number of molecules must be a concave function of the mole fraction X\ otherwise a 
double tangent construction is required to render the Gibbs free energy per molecule g concave. 24 * An 
interval for which ti g/dx 2 « 0 corresponds to phase separation. 
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corresponding value of the surfactant chemical potential /x s is equal to this minimum 
value; i.e. 

Ats = mm^s(M,N) = min(G^-MAto)/M (22) 

The proof is as follows. If /Lt s <niin (l s (M,N), then Glfj*>M(i 0 +Nii$ for all M 
and N. That is, the free energy of every aggregate exceeds the free energy of the 
monomers obtained by dispersing the aggregate. But if fi s > min iisiM^ N) 9 then 
G°m 9 s < Mfi 0 + Nfis'tor some M, N. The formation of aggregates with this M and N 
would lower the free energy. The problem is then reduced to finding min (is(M* JV). 

We shall see in the following sections how curvature energy and packing 
determine which aggregates (normal or inverse) will form. The formation of a 
microemulsion depends on the competition between micellar structures which 
absorb little oil and large "swollen" micelles which contain a large amount of oil. 

6.2 PACKING IN MICROEMULSION DROPS 

As in the case of surfactants we require now a model of microemulsion drops 
which will allow us to write down an explicit expression for /*s(Af, N). The simplest 
model of an aggregate is represented in fig. 4. For the time being we admit no specific 




Fig. 4. — Geometry of microemulsion drops; no penetration model. 



interaction between the oil and surfactant chains. The generalization of fi^ dis- 
cussed in section 2 is 

GUN = Mtio(brik)+^ (23) 

where a-A/N is the total area of the drop divided by the number of surfactant 
molecules, and ^o(bulk) is the chemical potential of oil (at the external pressure of 
the solution). It follows that the quantity to be minimized is 

Ms (A/,N) = j^ (24) 

where fi 0 is the chemical potential of the oil monomers. We restrict the argument to 
the case of excess oil for which fi 0 (b\xlk) == /i 0 . We emphasize agaiin that the 
particular form of curvature energy is immaterial for the argument which follows. 
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The variables a and R are specified through the relations 

. V-Mv«+No-*Rl, (25) 

so that the minimum fis{M y N) can be obtained by minimizing with respect to Rw 
and a in succession subject to the packing constraint / < / c . The minimum at fixed a is 
given by the smallest value of i? w allowed by packing. This corresponds to a drop 
with / = / c and no interpenetration of oil and surfactant for i)/a 0 / c > 3 or a spherical 
micelle for vf aok < h For water in oil microemulsions a similar argument leads to a 
form for fi$(M, N) identical to eqn (24) but with mo replaced by Mw and D replaced 
by -D. In this case the minimum corresponds to the largest i? w allowed by packing, 
/ = l C y and again no interpenetration. 

If curvature is assumed to be small, we can take a = a 0 and consider optimal drop 
parameters as a function of vfa 0 l c . We have 

y|«o-«w|-Afe;" l*o-*wHc, 4^=JVfl 0 , (26) 



whence 



3v R% +i?o-Rw^"-Rw i 2 



2 =l+x+x\ x=-zr (27) 



a 0 / c 3i?w ' Rvf 

with solution 



4 



3v 3 1 



, A ^- (28) 

Thus v/a 0 l c determines x and hence R 0 and i? w are determined if / c is known 
separately. Since R 0 /Rw<l, for normal drops v/a 0 l c < 1. The volume fraction of 
oil within the drop is 

For inverted drops the same equations hold, but now Ro-Ry/-U and the 
condition for existence is v/a 0 l c > 1. For inverted drops 

tfwater^ ^ =p. (30) 
y water ~ K surf * 

Volume fractions and other parameters of microemulsion drops are tabulated in 
table 4. Note that for oil-in-water systems geometry dictates that microemulsions, as 
opposed to swollen micelles, begin at v/a 0 l c 2> 0.8 (in the sense that <k>a>?). The 
problem of encapsulating a substantial oil fraction is critically sensitive to v/aoh> 
The observation that for inverted drops the aggregate of minimum free energy 
corresponds to the largest i? w allowed by packing has an interesting consequence for 
concentrated dispersions. For a given surfactant which forms inverse structures 
(v/a 0 l c > 1), uptake of water is greater on addition of oil. This can be illustrated by 
the simple model calculation for inverted micelles discussed below. Column 5 of 
table 4 lists maximum uptake of water by inverse micellar solutions as calculated 
from this model. These results are obtained as follows. Consider an inverse micellar 
system to be a dispersion Of equal-sized spheres of water in surfactant such that the 
head groups of all surfactant molecules reside at a water interface! For simplicity 
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Table 4.— (a) optimal parameters of oil-in-water microemulsions 

ASSUMING a =a„, NO PENETRATION OF OIL INTO TAIL REGION 



v/a 0 l c 


x — Ro/Rw 


«w//e 




0.5 


0.36 


1.58 


0.05 


0.6 


0.52 


2.10 


0.14 


0.7 


0.66 


2.96 


0.29 


0.8 


0.78 


4.64 


0.48 


0.85 


0;84 


6.31 


0.60 


0.9 


0.89 


9.65 


0.72 


0.95 


. 0.95 


19.70 


0.86 


1 


1 


00 


1 



Ko//c = *w/'c-l 



Table 4.— (6) optimal parameters of water-in-oil microemulsions 

(a = fl 0 , NO PENETRATION) 



v/aoU 


x — R 0 f R w 


R w /I e 


0w»t«r 


<f> v (inverted micelles) 


1 


1 


00 


1 


0.40 


1.1 


1.10 


10.32 


0.76 


0.34 


1.2 


i;i9 


5.31 


0.60 


0.29 


1.3 


1.27 


3.64 


0.48 


0.26 


1.4 


1.36 


2.80 


0.40 


0.22 


1.5 


1.36 


2.80 


0.40 


0.20 


1.6 


1.51 


1.95 


0.29 


0.18 


1.8 


1.59 


1.71 


0.25 


0.14 


2.0 


1.79 


1.26 


0.57 


0.12 



ignore head-group size. The volume fraction of water cf> w and radius R w of the water 
spheres is easily shown to be <£ w = J? w /(3tVa 0 +i?w). Let <*max be the maximum 
distance of any point of the surfactant region from the centre of the nearest water 
sphere. Clearly d mBX ^ # w + k For a given volume fraction <£ w , the minimum value 
of d max obtains for a f.c.c, lattice and can be shown to be 

Hence the largest volume fraction allowed is that corresponding to a f .ex. lattice with 
rfmax + R w + /c Eliminating /? w » the tabulated results follow. 



6.3 EFFECT OF CURVATURE IN MICROEMULSION FORMATION (H.L.B.) 

If v f a 0 and l c can be estimated, table 4 gives some idea of the size of drops 
expected and therefore of design characteristics. The model is clearly deficient in at 
least two important respects. One is its lack of oil specificity which problem we defer. 
The other is the effect of curvature which can affect allowed optimal structures in 
some circumstances. A mathematical description of curvature effects which spans 
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the full gamut of expectations and trends at a qualitative level only can be formulated 
as follows: 

(1) For non-ionic surfactants we have for the surface contribution to pts(Af, N) 

d/rV 



^i a + a(l+°D/R))> (31) 



where, as for eqn (18H20) the bilayer optimal area can be characterized for 
convenience in terms of a capacitor model for head group interactions, ao- 
<J(27rq 2 D/sy). As for our analysis of vesicles, for definiteness only we take 
a 0 = V(372D/A 2 ). 

(2) An identical description can be used for ionic surfactants. Debye-Huckel 
theory leads again to the capacitor model with D = #c~ l , the Debye length, with 

~ekT n ° 

for a 1:1 electrolyte, n 0 is the bulk electrolyte concentration. Taking 5=78 
for water, q =4.8xl(T 10 e.s.u., fer=4.116xl0~ l4 erg at 25°C, we have 
£> = 3.03<r 1/2 L where c is the salt concentration in moldm"" 3 . Taking further 
y = 50 dyn cm" 1 as the surface tension of the hydrocarbon-water interface for rough 
estimates we have the correlation between salt concentration and head-group area 
indicated in table 5. The import of table 5 is that for ionic surfactants, salt 
concentration or composition must provide a substantial contribution to head-group 
curvature and can markedly affect conditions under which microemulsions can form. 



Table 5.— Effect of salt concentration on bilayer head-group area 



c/mol dm -3 


D/A 




»/a 0 /c(»// c =30A 2 ) 


1 


3.03 


33 


0.893 


0.5 


4.29 


40 


0.751 


0.1 


10.0 


60 


0.502 


0.05 


13.6 


71 


0.422 


0.01 


30.3 


106 


0.282 


0.001 


96 


189 


0.159 



Clearly both models of curvature encompassed in eqn (27) are extreme and are 
not to be taken at a quantitative level. However, characterization of curvature 
through such models may be useful in two respects: (a) Assignment of D to a 
particular non-ionic surf actant/water system can be used to predict behaviour for 
others in a homologous series (e.g. C„EO m ); similarly for ionic surfactants with 
varying salt concentrations; (6) The very extremity of such models is useful if, given 
such flexibility, the substantial conclusions of the z&roth-order theory are confirmed. 

We consider then the form eqn (31) and optimize /i% with respect to a, taking 
R = 3v/a for spherical micelles and fl w = R - / c /[|- V(3u/a/ c -!)] [eqn (28)] for 
microemulsion drops. Differentiating eqn (31) we have 

* al R /„ a 3R\ 



D. J. MITCHELL AND B. W. NINHAM 619 

where 

.—■—=-1 spherical micelles (R - 3v/d) (33a) 

a BR 3(v/al e ){R/Q . . • . /<Mt < 

jri.w- ,(,, ™ ,,d ^ (33 * } 

For v/al c > 3 we have drops and for v/al e < 3 micelles. At v/al e - 3, R = / e , eqn (33ft) 
reduces to 

1? 3a ~ 

so that dfi%/da is continuous. Thus we can find the optimum structure if it exists by 
finding the zero of 8n%/da. The same equations hold for inverted drops if R is 
replaced by -R. The equations have been solved numerically for a range of 
surfactant parameters. Typical results are displayed in table 6. 

Table 6.— Effects of head-group curvature on parameters of 

microemulsion drops 
(i) k - 17.5 A, V = 1060 A 3 (typical double-chained surfactant) 



D/A 


ajA 2 


a/A 2 


*w/A 


Ri/A 


v/a 0 l t 


v/al c 


1 


19.3 


21 


-13.5 


-13 


3.1 


2.9 


3 


33 4 


38.5 


-36.5 


-27 


1.8 


1.56 


5 


43 


51 


-102 


-50 


1.4 


1.2 


7 


51 


60 


+3663 


-105 


1.18 


0.96 


9 


58 


68 


144 


-479 


1.04 


0.86 


11 


64 


74 


85 


276 


0.94 


0.81 


13 


69.5 


79.5 


65 


121 


0.86 


0.75 


17 


80 


89 


47 


62 


0.74 


0.67 


21 


88 


94 


42 


48 


0.68 


0.64 


30 


106 


104 


34.5 


33 


0.57 


0.57 


43 


126 


112. 


30.5 


26 


0.47 


.0.54 


96 


189 


124 


26 


16 


0.32 


0.48 


303 


336 


132 


24 


9 


0.18 


0.45 


(H)/c = 


17.5^0 


= 525 A 3 (TYPICAL SINGLE-CHAIN SURFACTANT) 


D/A 


ao/A 2 


a/A 2 


i?w/A 


R°JA 


v/a 0 tc 


v/al e 


3 


33 


35 


107 


165 


0.90 


0.84 


9 


56 


58 


28.5 


29 


0.52 


0.51 


13 


69.5 


67 


.24 


23 


0.43 


0.45 


17 


81 


73 


21 


19 


0.37 


0.41 


21 


88 


77 


20.5 


18 


0.34 


0.39 


30 


105 


84 


. 19 


15 


0.28 


0.36 


43. 


126 


89 


18 


12 


0.24 


0.34 


76 


189 


97 


16 


8 


0.16 


0.30 


303 


336 


103 


15 


4.7 


0.09 


0.29 



Comparison of raicroemulsion drop parameters as given by zeroth-order theory (a = a 0 ) 
with predictions of capacitor model for head-group curvature. Negative radii correspond to 
inverted drops. R w denotes R w as computed by zeroth-order theory. 



620 



MICELLES, VESICLES AND MICROEMULSIONS 



Several remarks are in order. Recall that if curvature is ignored, the zeroth-order 
approximation predicts (i) micelles v/a 0 l G <h GO normal drops 3<v/a 0 l c <l; (Hi) 
inverted drops v/a 0 l c > 1 . By comparison, when curvature is taken into account, the 
situation is modified. Now we predict normal drops in a region where formerly we 
would have expected inverted drops. For smaller values of v/a 0 l c curvature effects 
give larger drops, indeed can lead to drops where zeroth-order theory gives micelles- 
Note that for single-chained surfactants inverted drops only occur for low values of 
ao, i.e. very weak head-group interactions, The indications are that curvature plays a 
larger role in the determination of the dimensions of microemulsion drops than for 
vesicles, where there is a partial cancellation of curvature effects. 

It can easily be shown from eqn (32)-(33) that the condition for existence of 
inverted drops is 

«/a 0 f c >V(l+D/Q, 

and for micelles 

i jq+2P/m 

0/ao/c< 3 (l+2>/4) • . 

Note that the zeroth-order theory is in rough qualitative agreement. In essence 
then, these calculations quantify the older notion of hydrophilic-lipophilic balance 
(h.l.b.). We shall expand on this in the following section. 

7. DISCUSSION AND CONCLUSION 

7.1 RECAPITULATION 

We are now in a position to put the whole story together. To recapitulate: We 
have presented a simple model of surfactant aggregation in which the aggregate size, 
shape and composition is determined essentially from the packing ratio t>/a 0 f c - 
Components of this apparently magic number are v, the partial molecular volume of 
surfactant, a 0 , the head-group area of a surfactant molecule and / c , which is the 
maximum chain length. The head-group area a 0 is assumed to be determined by the 
competition between head-group interactions and the interfacial tension between 
the hydrocarbon tails and water. It is thus primarily a property of the head group. 
The more "hydrophilic" the head group the larger a 0 will be. The parameter l c is 
assumed to be a property primarily of the tail. It is taken to be the maximum length 
of the tail in the sense that a further increase in its length will result in a large increase 
in free energy. 

We have generally assumed that the tails are completely flexible and liquid-like 
provided their length / is I$ss than / c , and that / is constrained to be less than / c . The 
aggregate then is assumed to adopt the optimum head-group area. Then the 
head-group curvature energy, which is treated as a perturbation, is assumed to shrink 
the aggregate to the smallest aggregate satisfying the above constraints. The smallest 
aggregate, for v/a 0 l c >% will have / * / c . However, the situation would not change 
significantly (except for vesicles) if 7 C represented an optimum length rather than a 
maximum length. It is essentially the tail length at which the tail free energy resists 
further lengthening under the influence of head-group curvature energy. 

In the case of vesicles we have found it necessary to introduce a lower limit al c on 
the tail length. The free energy of the tails is assumed to rise rapidly if 7 falls below 
al c . However, it would appear that tail lengths can vary by at least a factor of two. 19 
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The main omission in this lowest level theory is the possible effect of curvature on 
the free energy of both head group and tail. We have explored this matter in sections 
4.2, 5.2 and 6.2. 



7;2 H.L.B., AND EFFECTS OF SALT, TEMPERATURE. AND 
COSURFACTANTS 

The packing ratio v/a Q l c provides a measure of the hydrophilio-Iipophilic balance 
(h.l.b.). For values of v/a 0 l c < 1 we predict "normal", i.e. convex, aggregates and for 
v/a 0 l c > l we predict inverse micelles and drops. That is, for i>/a 0 / c < 1 we predict a 
tendency to curve spontaneously into normal structures and for v/aoU> 1 to curve 
spontaneously into inverse structures. 

The packing ratio is affected by many factors including hydrophilicity of head 
group, ionic strength of solution, pH, temperature and the addition of lipophilic 
compounds such as cosurfactants. The more hydrophilic the head group the stronger 
the repulsion between the head group and the larger a 0 will be. The smaller the value 
for v/a 0 h> the smaller will be the micelles formed. Thus (single chained) ionic 
surfactants (at low salt concentrations) tend to form small spherical or globular 
micelles 1 whereas monoglycerides tend not to form micelles at all but form vesicles or 
liposomes or lamellar particles in dilute solution. 12,16 The non-ionic surfactants 
polyoxyethylene dodecyl ether (Ci 2 E m ) form micelles at room temperature if 
their head groups are long enough (e.g. Ci 2 E 6 , Ci 2 E 8 ) but C i2 E 4 forms lamellar 
dispersions. 12 

The addition of NaCl to micellar solutions of sodium dodecylsulphate causes the 
micelles to grow in size and become long rod-shaped micelles. 1,28 Similarly the 
addition of HC1 to micellar solutions of potassium oleate causes the micelles to grow, 
and, if the pH is lowered sufficiently, the solution goes cloudy and contains 
liposomes. 25 In both cases the repulsion between head groups decreases, Le. a Q 
decreases and v/aoh increases, so that our theory would predict a change from 
spherical micelles to cylindrical micelles and eventually to vesicles. 

The incorporation of cosurfactants into surfactant aggregates would be expected 
to increase the mean volume per surfactant molecule without affecting appreciably 
either a 0 or / c . 26 Consequently v(a<J c would increase with addition of cosurfactant, 
leading successively to an increase in micelle size, the formation of long rod-shaped 
micelles and eventually the formation of lamellar phase or lamellar dispersion. 26 ' 27 
The theory also predicts the observed increase in stabilization of hydrocarbons by 
ionic surfactant solutions due to the addition of cosurfactants (e.g. cf. Shinoda 22 ). 
There is no need to invoke specific interactions (which are discounted by Shinoda and 
Friberg, 21 and by Wennerstrom and Lindman 26,28 ). 

The effect of temperature on v/aok is difficult to predict without a better 
understanding of the forces at play. However, many of the observed phenomena for 
biological lipids 2 - 12 - 21 as well as 1221 for non-ionic surfactants such as polyoxy- 
ethylene) alkyl ethers may be accounted for if we assume v/a 0 U increases with 
temperature. Micellar size would be expected to increase with temperature. Non- 
ionic surfactants do tend to go cloudy when the temperature is raised above the aptly 
named cloudy point 22 indicating the formation of large aggregates [or flocculation of 
aggregates, the phase diagrams of poly(oxyethylene) alkyl ethers display a lower 
consulate curve 22 near the cloud point which indicates an attractive interaction 
between the aggregates above the critical point]. The surfactant Ci 2 E4 forms 12 a 
lamellar dispersion at still higher temperatures. 
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If vjaoh increases with temperature one would expect the solubilisation of 
hydrocarbons in non-ionic surfactants to increase with temperature (as observed) 22 
until u/Wc reached the value 1 where phase inversion would be expected. At higher 
temperatures water-in-oil microemulsions would be expected and the solubilisation 
of water would decrease as temperature rises (again as observed). 21 ' 22 

As ao increases with increase in length of the poly(oxyethylene) head group we 
would expect then an increase in the cloud point and phase-inversion temperature 
(p.i.t.) with head-group length (as observed)/ 2 

73 OIL SPECIFICITY 

The model of microemulsion drops presented in section 6 ignored the possibility 
that the oil might penetrate into the surfactant interphase. Some phenomena of oil 
solubilisation appear to require the assumption that hydrocarbons, some more than 
others, do penetrate into the interphase. As a consequence, the packing ratio v/a 0 l c 
is effectively increased, due to an increase in the interphase volume per surfactant 
molecule. 

Experiments on the solubility of hydrocarbons in lipid lamellar phases 29,30 
suggest that smaller alkanes penetrate into the surfactant layers more than larger 
alkanes. If we assume this to be true then this accounts for the variation with 
alkane chain length of cloud points of non-ionic surfactant solutions with saturated 
hydrocarbon. 22 This observation can be illustrated and quantified through the 
following example. From the experiments discussed by Gruen and Haydon 29 we see 
that decane swells the hydrocarbon layer of a glycerol monoleate black film to double 
its thickness (from 24 to 48 A); Experimentally it is known that the head-group area 
remains constant. Thus the volume per surfactant molecule doubles and so v/a 0 l c 
also doubles. On the other hand, the p.i.t. (which corresponds to i>/ao/ c = 1) of the 
decane-Ci 2 E4-water system 31 is approximately 20 °C which is coincidentally the 
onset temperature of lamellar dispersion for aqueous solutions of C 12 E 4 12 cor- 
responding to v/a 0 l c =-h Thus also in this example decane appears to swell the 
surfactant hydrocarbon tail volume by a factor of two. With longer chain hydro- 
carbons v/a 0 l c changes less markedly 29 - 30 (e.g. for hexadecane by a factor of f)» 

Long-chain alkanes elevate the cloud point. 22 This is to be expected since v/a 0 l c 
is hardly altered, and our model predicts that large, non-spherical micelles will be 
replaced by small microemulsion drops. 

7.4 THE INTERFACIAL TENSION CONTROVERSY AND 
THE SURFACTANT PHASE 

Some authors 32 have suggested that a microemulsion forms spontaneously as a 
result of the interfacial tension between the oil and water phases becoming zero or 
even taking negative values. If, on addition of surfactant or cosurfactant to the 
system, the interfacial tension becomes negative the system would tend to increase 
the area of the oil-water interface, taking up surfactant from the solution, until the 
interfacial tension returned to zero. The system would then be stable because: If the 
area were to decrease, the interfacial tension would become negative, favouring an 
increase in the area, whereas if the area were to increase the interfacial tension would 
become positive favouring a decrease in area. However, if we ignore such factors as 
curvature energy, the system would be indifferent to whether an oil-in-water or 
water-in-oil emulsion formed, or to what size distribution of drops formed, provided 
only that the total area achieved the required value. 
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By contrast, if we invoke a curvature energy to determine the size of the drops 
and to determine which phase is dispersed in the other, we find that the micro- 
emulsion will form while the tension of the oil-water interface is still slightly positive. 

This can be seen from the following argument: If curvature energy favours a drop 
of finite, size it will do so only because this drop has the minimum value of 
(GW.n -Mti 0 )/N or of (Gm.n .— Mfiw)/N. The drops will form when the chemical 
potential of the surfactant is equal to this minimum value. On the other hand, for an 
infinite drop of oil {Gm.n -Mfio)/N represents the interfacial Gibbs free energy per 
surfactant molecule and if the interfacial tension of the infinite drop is zero, this 
interfacial Gibbs free energy is also equal to the chemical potential of the surfactant. 
Thus the finite drop will form at a surfactant chemical potential lower than the 
chemical potential at zero interfacial tension, i.e. when the interfacial tension is 
positive. 

Shinoda and Friberg 21 sketch the variation of interfacial tension between the oil 
and water phases with temperature for a system containing equal amounts of water 
and oil [47.5 wt %, and 5 wt % of poly(oxyethylene)nonylphenyl ether]. The inter- 
facial tension drops to a very low value at the p.i.t. This might be understood roughly 
from our model if we assume again that v/a 0 l c increases with temperature, taking on 
the value unity at the p.i.t. Below the p.i.t. v/a 0 l c < 1, curvature energy favours an oil 
water dispersion and this implies a slightly positive interfacial tension. At the p.i.t. 
v/a 0 l 9 =l, curvature energy favours zero curvature leading to a zero interfacial 
tension. At and near the p.i.t. therefore one expects either large drops with little 
curvature or an essentially lamellar structure, or some such bicontinuous phase. 37,38 
Such large structures would be expected to interact strongly and form a separate 
phase; hence the surfactant phase which is commonly observed 21 at the p.i.t. The 
above argument and the experimental data suggest then that not only can micro- 
emulsions form at positive interfacial tension, but that a positive tension is in fact 
required. Very small interfacial tensions will occur only when the surfactant phase 
forms. 

7.5 INTERACTIONS 

So far we have not taken interactions between aggregates into account. The 
theory has been restricted to dilute oil-in-water dispersions in the presence of excess 
oil, or vice versa. We have correlated the type of microemulsion which forms to the 
type of aggregate which forms in the absence of oil. If we ignore the penetration of 
oil into the surfactant tails we would predict that surfactants forming normal micelles 
or vesicles should, on addition of oil, form small swollen micelles or microemulsion 
drops; and that surfactants which form inverse micelles should form water-in-oil 
dispersions on the addition of oil. The type of aggregates formed by surfactant 
solutions can change from normal to inverse structures on increase in temperature or 
addition of cosurfactant. Penetration of oil into the surfactant tails will shift the 
balance towards inverse structures. 

Changing the water content of a surfactant solution can also change the type of 
structure which is observed. 12,27 There are a number of possible related reasons for 
this. . Firstly, for a fixed volume fraction the distance between aggregates can be 
increased by decreasing the area per head group which, through the packing 
conditions, can change the shape; from sphere to cylinder or from cylinder to 
lamellar. Increasing the distance between aggregates may decrease the (repulsive 
hydration or ionic) interaction energy. 33 Secondly, the distance between aggregates 
can be increased by changing the shape from spherical to cylindrical or from 
cylindrical to lamellar, thus decreasing the interaction energy at the expense of 
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curvature energy. Thirdly, decrease in water content may reduce the amount of 
water available for hydration of the head groups. This would lead to a smaller 
head-group area, larger value of v/a Q l c and thereby could possibly induce a change in 
shape of the structures present Whatever the reasons, change in structure is 
observed in surfactant solutions with change in the concentration. One might 
reasonably expect that analogously to the dilute solutions considered in this paper, 
the type of microemulsion formed on the addition of oil might be correlated with the 
structures present in the absence of oil, taking into account the effects of penetration 
of oil into the surfactant tails. Since surfactant solutions can change from normal 
micellar solutions to lamellar phase to inverse micellar solutions, on decreasing the 
water content, one might reasonably expect that water-oil-surf actant-cosurfactant 
systems may change from oil-water microemulsions to water-oil microemulsions 
via a lamellar phase as is indeed commonly observed. 34 Some surfactant solutions 
form reversed hexagonal phases 27 so that it is not surprising to find reverse hexagonal 
phases in oil-water-surf actant-cosurfactant systems. 34 

We are well aware of controversies and debate on opposing points of view. 
Although the older mixed-film theory has merit and embraces part of the story, 
Friberg and Buraczenska 34 have demonstrated that the distinction between micro- 
emulsions and swollen micelles is artificial. Our object in this paper has been to 
develop a unified approach to the subject of microemulsions which also embraces 
micellar systems. On that unremarkable note we rest our case. 

We owe a debt to several people: The problem was posed to us by J. Th. 
G. Overbeek. Extensive discussions with G. J. T. Tiddy, our colleagues J. N. Israel- 
achvili, S. Marcelja and B. A. Pailthorpe are gratefully acknowledged. We are 
also much indebted to S. Friberg and K. Shinoda not just for the care and thorough- 
ness of their experiments, but also because of the clarity of their expositions on the 
subject. B.W.N, is grateful to S. Friberg for much encouragement. 

APPENDIX 

THE LAW OF MASS ACTION AND THE CONCEPT OF AN AGGREGATE 

As already emphasized in section 2, the whole subject of association colloids is 
based on intuitive notions and is fraught with difficulties which devolve on 
definitions. The treatment in this paper is largely centred around the law of mass 
action and relies heavily on the notion of an aggregate (micelle, vesicle, drop, etc.). 
Here we wish to explore the difficulties and limitations involved in these concepts and 
the nature of the approximations made. 

The concept of an aggregate is not precisely defined. There is a certain arbitrari- 
ness about any mathematically precise definition but any definition adopted must 
conform closely with our intuitive concepts and the predictions of the theory must 
not depend critically on our (arbitrary) choice of definition. 

(1) a one-dimensional illustration. The difficulties alluded to are perhaps 
most clearly illustrated by means of a simple one-dimensional model which allows 
analytic solution. Consider a one-dimensional gas with nearest-neighbour inter- 
actions. Its solution 35 is 




'o 



(Al) 
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where /x is the chemical potential, /3 = l/kT> p is the pressure and u(x) is the 
potential of interaction between nearest neighbours. [We assume that the distance 
x is normalized by some suitable length such as the diameter of molecules so that 
the integral in eqn (Al) is dimensionless. The pressure p is then also scaled 
appropriately.] Suppose that 

If* 

u(x) = Q iorx>R and — expH?«(*)]dx»l. (A2) 

XV Jo 



Then for QpR « 1, eqn (Al) can be approximated by 

• R 

exp ' 



(-fti)-f c«p[-^(x)]dx+i . (A3) 
Jo (3p 



whence 



where 



I=\ exp [-/?«(*)] dx. (A5) 
Jo 



>R 

Thus the concentration p given through the thermodynamic relation p -dp/dfi is 

p= £ m^mplfiNM. (A6) 

By comparison the law of mass action, eqn (1), gives 

P=I^=I NtxpQ3NiJL)exp(-pNn$). (A7) 
jv-i N-i 

Eqn (Al) and (A7) are equivalent if we make the identification 

exp(-^) = /^ 1 =J o .dx 1 ^ dx 2 .. djc^iexp[^ !«(*/)]. (A8) 

We have shown then that, under the conditions 

iSpJR « 1, -J o exp[-^U)]djc»l, (A9) 

the exact thermodynamic equation (Al) reduces approximately to the law of mass 
action with p# determined by eqn (AS). In eqn (A8) exp is seen to be the 

partition function of an AT-particle system forming an aggregate where the particles 
are constrained so that no two nearest neighbours are more than a distance R apart. 
The value of R is arbitrary so long as it satisfies condition (A2), but the conclusion is 
insensitive to the choice of R. The condition for the above concept of aggregate to be 
meaningful is that the conditions (A2) and (A8) should hold for some R> so that the 
potential u (x) should be strongly attractive and short-ranged. 

(2) three-dimensional, one-component; For a three-dimensional system the 
justification of the law of mass action is rather more complicated. Consider first a 
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one-component system, Le. ignore any solvent The thermodynamic potential 
fl= -PV is related to the grand partition function by 

exp (fiPV) = Sfo, V) = 1 exp (fifiN)Z N (V) f (A10) 

where V) is the canonical partition function which is defined by 

Z(V) j^d* . „ J dr* exp HMn, . . . , #*)] (All) 

where A = (2irh 2 /mkT) m is the thermal wavelength, r u r 2% . . . , r# are the coor- 
dinates of the N particles and the Hamiitonian is taken to be 

H„=Z PV2m + U N (r u r 2t ... f rs). (A12) 

Assume now that some arbitrary criterion has been chosen for defining an 
aggregate* For example we might agree that m molecules form an aggregate of 
aggregation number m (an m-mer) provided that no molecule of the aggregate is 
more than some prescribed distance R from some other molecule of the aggregate 
and no subgroup of molecules is further than R from the rest of the aggregate, and 
provided of course that the m molecules do not belong to a larger aggregate. Then 
for each configuration (r u r 29 . . . , the N molecules partition into aggregates. Let 
n m denote the number of m-mers for m 555 1, 2, 3, . . . , AT. The partition function 
Zn( V) can then be rewritten as a sum over all configurations consistent with the 
partition, or 

^"An^^ nn inrmn"mi f d 3i ^ exp [-^ (r u . • . , r N )], (A13) 

where $l{{n m }) denotes the region of configuration space consistent with the parti- 
tion into rti monomers, n 2 dimers, . . . n m m-mers, etc. We have included in eqn 
(A13) a combinatorial factor which gives the number of ways of partitioning N 
indistinguishable molecules into n m aggregates of size m, m = 1, 2, . . . , N. Any two 
configurations which differ only by the interchange of two molecules should be 
considered identical. The necessary approximations consist now in ignoring inter- 
actions between aggregates and extending the range of integration in eqn (A12) 
to allow aggregates to move through each other. The errors induced by these 
approximations is negligible, and allows a desirable simplification. Thus we have 

^W-j^ m /(Ai4) 

or 

^(v)4in§)"" ,(Ai5) 

where 

/« = drx...dr^Qxpt-0u m (r u ...,r m )] • (A16) 
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the integration being oyer all configurations satisfying the constraints imposed by the 
definition of an aggregate. If we recall that the set n m are. subject to the constraint 
X m mn m = N, we have from eqn (A10) 

m » m =o n m \ \ A / {ml) 

UVjtp^y^V). (AH) 

Thus 

p\ 3 =ZTnexp(pmri(—k^): ' (A18) 

This is equivalent to the law of mass action provided we take 

exp(^m^ m ) = — yf-j. (A19) 

For the above result to have any validity it is necessary not only that the system be 
sufficiently dilute that interactions between aggregates (and monomers) be negligible 
but also that I m be independent of the precise definition of an aggregate. This 
requires that the potential of interaction be short-ranged and strongly attractive. 
The thermal wavelength A does not play any fundamental role in the above argument 
except (primarily) as a length scale to keep the above equations dimensionally 
correct. Any other suitable scaling length might be chosen, e.g. the diameter of a 
molecule (cr). If we rewrite eqn (A18) and (A19) 

pcr 3 ^Zexp08m^)( w ^-3) (A20) 
exp {-pmtxli) - . 3m-3 , (A21) 

ttt\(T 

we will have merely replaced fi by ft + kT In (A/<r) 3 and mfin by m/t^-f 
(m \)kT In (A/tr) 3 . That is we have merely changed concentration units and 
changed the reference state with respect to which the chemical potential is defined. 
Thus the classical potential of the ideal gas is now kT In (per 3 ) instead of kT In (pA 3 ). 
The reference state is the hypothetical ideal gas of density p = 1/cr 3 . This change is 
convenient wheii working at a fixed temperature but is awkward when considering 
changes in temperature since the reference state is temperature dependent. 

(3) three-dimensional, MULTTCOMPONENT systems. The above argument can 
be extended to a multicomponeht system. For example consider an aggregating 
solute B dissolved in a solvent A. We may write the grand partition function as 

(A22) 

where 0 = 1 /fcT, P is the pressure, V the volume, ja a and (i B the chemical potentials 
of solvent and solute, respectively, A a and A B the thermal wavelengths and 
Us At N B (^r) the potential energy of the N A solvent molecules and N B solute 
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molecules whose positions are given by q and r, respectively. If we let P 0 denote the 
pressure of pure solvent at the given (ia and T then 

exp (fiPoV) = I ^H7$r f d '" A ex P H^ofo r)l (A23) 
We define the excess free energy u NB (r) of the N s solute molecules as 

exp[-iSfi NB (r)] = exp(-^PoV) £ ^TT^ f d ^ exp r)1 

. (A24) 

so that eqn (A22) may be rewritten as 

exp (faV)=Z eX * ( f *% b) [ d 3 % exp [^ B (r)], (A25) 

where tt =P-jP 0 is the osmotic pressure of the solution relative to pure solvent at 
pressure P and temperature T. Eqn (A25) is completely analogous to eqn (A10) and 
( Al 1) with the gas pressure P replaced by the osmotic pressure n and the potential 
energy of the gas molecules replaced by the excess free energy of the solute 
molecules. The same argument may therefore be applied to obtain the law of mass 
action for the solute and obtain formal expressions for the /iff- The /t ^ are thus given 
as functions of fi A and T but for dilute solutions (i A is approximately that of the pure 
solvent at the appropriate pressure (usually 1 atm). 

So far we have only obtained a formal expression for (il?. In the next subsection 
we shall attempt to obtain a physical interpretation for pfc 

(4) physical meaning of fis. We have obtained a formal expression [eqn (A19)] 
for Mm in terms of the partition function I m defined by eqn (A16). The partition 
function I m resembles a canonical partition function except that for our integral I m is 
over configurations which constitute an aggregate whereas the canonical partition 
function is an integral over configurations in which the molecules are confined to a 
fixed volume. However, as is well known in statistical mechanics, the major 
contribution to the canonical partition function comes from configurations with a 
uniform density filling the whole container. One might likewise expect that the 
major contribution to J m comes from configurations with a uniform density equal to 
the average density of an aggregate and that I m would be approximately equal to the 
canonical partition function of a system of m particles with this average density. The 
mfi ° m would be approximately the Helmholtz free energy of a system of m molecules 
of this average density. Since the system is of finite size it will have in addition to a 
bulk free energy term a surface energy and possibly curvature energy. 

It should be clear from the above example that the constraints defining an 
aggregate should confine the molecules to a sufficiently small volume so that a natural : 
aggregate will form spontaneously just as a gas if confined to a small volume will 
collapse to form a liquid drop. 

The interpretation of fi° m for a solution is more involved. The expression [eqn 
(A 16)] for I m will involve the excess free energy defined by eqn (A24). If we assume 
phase separation of the solute and the solvent then clearly mti° m equals the bulk 
Helmholtz free energy for m molecules at the average aggregate density +P 0 V &ug + 
surface terms + curvature terms. If the aggregate is incompressible the sum of the 
first two terms is the bulk Gibbs free energy of m molecules at pressure P 0 [i.e. at the 
external pressure (usually 1 atm) not at the internal pressure of the aggregate which 



D. J. MITCHELL AND B. W. NINHAM 629 



may differ from P 0 by the Laplace pressure 2y/r where y is the interf acial tension and 
r the radius of curvature of the drop], 
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